Astronomy & Astrophysics manuscript no. p6v3 


© ESO 2009 


May 21, 2009 





Radiative transfer in cylindrical threads with incident radiation 

VI. A hydrogen plus helium system 

P. Gouttebroze 1 and N. Labrosse 2 

1 Institut d'Astrophysique Spatiale, Univ. Paris XI/CNRS, Bat. 121, F-91405 Orsay cedex, France 
e-mail: pierre .gouttebrozeOias .u-psud. f r 

2 Department of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, Scotland 
e-mail: n. labrosseOphysics . gla. ac .uk 

Received 8 December 2008 / Accepted 

ABSTRACT 

Context. Spectral lines of helium are commonly observed on the Sun. These observations contain important informations 
about physical conditions and He/H abundance variations within solar outer structures. 

Aims. The modeling of chromospheric and coronal loop-like structures visible in hydrogen and helium lines requires the 
use of appropriate diagnostic tools based on NLTE radiative tranfer in cylindrical geometry. 

Methods. We use iterative numerical methods to solve the equations of NLTE radiative transfer and statistical equilib- 
rium of atomic level populations. These equations are solved alternatively for the hydrogen and helium atoms, using 
cylindrical coordinates and prescribed solar incident radiation. Electron density is determined by the ionization equi- 
libria of both atoms. Two-dimension effects are included. 

Results. The mechanisms of formation of the principal helium lines are analyzed and the sources of emission inside the 
cylinder are located. The variations of spectral line intensities with temperature, pressure, and helium abundance, are 
studied. 

Conclusions. The simultaneous computation of hydrogen and helium lines, performed by the new numerical code, allows 
the construction of loop models including an extended range of temperatures. 

Key words. Methods: numerical - Radiative transfer - Line: formation - Line: profiles - Sun: chromosphere - Sun: 
corona 



1. Introduction 

Observation of the upper solar atmosphere with high an- 
gular resolution reveals a wealth of filamentary structures 
produced by magnetic fields. For modeling these objects, we 
developed a series of NLTE radiative transfer codes that are 
described in the present series of papers. Among the fila- 
mentary objects relevant to this kind of modeling, we can 
mention: 

— cool coronal loops, 

— chromospheric fine structure (cf. Patsourakos et al. 
[2007]), 

— prominence (or filament) threads (cf. Heinzel [2007J) , 

— spicules. 

Paper 1 (Gouttebroze I2004|) was dealing with ID (i. e. 
radius dependent) cylindrical models, a case applicable 
to cylindrical structures with vertical axis, exposed to an 
incident radiation field independent of azimuth. Papers II 
and III (Gouttebroze 12005] and |2006J respectively) treated 
the case of 2D (radius and azimuth dependent) cylinders. 
Paper II was restricted to a 2-level atom, while Paper III 
used a multilevel hydrogen atom. Papers IV and V 
(Gouttebroze 120071 and I2008[) were dedicated to radiative 
equilibrium and velocity fields, respectively. All these 
papers were dealing with the hydrogen atom. Certain 
amount of helium was included in the state equation, but 
it was assumed to be neutral and without any influence on 



the radiation field, as well as on the electron density. 

In the present paper, we assume that the cylinders 
are filled with a mixture of hydrogen and helium, and 
treat NLTE radiative transfer and statistical equilibrium 
of level populations for both atoms in two dimensions. The 
helium model atom includes the three stages of ionization, 
and the electron density is recomputed at each iteration 
in order to satisfy the equation of electric neutrality. In 
Sect. 2, we describe the computational methods used in 
the new numerical code. The results concerning hydrogen 
and helium ionization are detailed in Sect. 3. In Sect. 4, 
we study the formation of helium lines using a reference 
model defined in Sect. 2. Finally, in Sect. 5, we show how 
the helium line intensities react to changes in temperature, 
pressure and helium abundance. 

2. Numerical methods 

2.1. Formulation 

The computation includes the numerical solution of the 
equations of NLTE radiative transfer for hydrogen and he- 
lium atoms, statistical equilibrium of level populations (for 
both atoms), pressure equilibrium, and electric neutrality. 
As in Papers II and III, the object under consideration is a 
cylinder of diameter D, whose axis makes an angle a with 
the local vertical to the solar surface (this angle may vary 
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between and 90°). The two active dimensions for radia- 
tive transfer are the distance to axis r and the azimuth 
i/j. The method of resolution is of the MALI type (Rybicki 
and Hummer 1991). The special form of equations for two- 
dimensional azimuth dependent (2DAD) cylindrical geom- 
etry and the method of solution, which are described in 
details in Paper II, will not be repeated here. The equa- 
tions of statistical equilibrium, independent of geometry, 
are treated in Paper I. In Paper II, we also described the 
method to compute the intensities incident on the cylin- 
der, at different wavelengths, from the knowledge of the 
emission by the Sun, the inclination a, and the altitude H . 
These incident intensities are also functions of azimuth, ex- 
cept in the special case a = 0. The condition for electric 
neutrality may be written 



A c = Nmi + Nueii + 2JV, 



VHcIII 



(1) 



where A e is the electron density, -/Vhii the number density 
of ionized hydrogen (protons), and iVnoii an d Ahciii the 
number densities of helium atoms in the second and third 
stages of ionization, respectively. If Ah and An e are the 
total densities of hydrogen and helium, respectively, and 
JVhi and An e i the corresponding densities of neutral atoms, 
the law of particle conservation yields 



A H = A H i + Ann 
and 



N 



He 



iVl 



Hcl 



A 



Hell 



A 



HcIII 



Besides, the gas pressure is 



P g = (A H + A Ho + A ) k T , 



(2) 



(3) 



(4) 



where k is the Boltzmann constant and T the temperature. 
The parameter to check the convergence of iterations is 
the electron-to- hydrogen ratio lo = (A /Ah). If An c is the 
(He/H) abundance ratio, Eq. (4) becomes 



P g = A H (l + ^Hc + ^) kT 



(5) 



The model atom for hydrogen includes 5 discrete levels plus 
1 continuum. The equations of radiative transfer for the 
10 discrete transitions and the 5 bound-free transitions are 
treated in details, except in the case where the optical thick- 
nesses are very low, which generally happens for continua 
from subordinate levels. The model atom for helium is the 
same as in Labrosse and Gouttebroze (2001, 2004]). It con- 
tains 34 levels: 29 for He I, 4 for He II, and 1 for He III. 
The main part of this model comes from the neutral he- 
lium model of Benjamin et al. (1999) . It is complemented 
with parameters of various origins, as described in Labrosse 
and Gouttebroze (|200l|) . The number of permitted radia- 
tive transitions is 76, but most of them arc optically thin 
under usual conditions. All discrete transitions, for hydro- 
gen as for helium, are treated under the assumption of com- 
plete frequency redistribution. 

2.2. Computational scheme 

The computation is organized along two parallel series of 
routines, the one for hydrogen, the other for helium. In 
each series of routines, the variables dependent on the 
atomic structures as atomic parameters, populations, in- 
tensities in different transitions, are gathered into a spe- 
cific COMMON, which is ignored by the main program. 



This main program only contains geometrical (D, r, tp, etc.) 
and physical (P g , T, vt, etc.) variables and the populations 
(Ahii, AhoII, etc.) necessary to determine the ionization. 
The main phases of computations are: 

— Initialization: determination of geometrical, physical 
and atomic parameters. The incident intensities are also 
computed for each position at the surface of the cylinder 
and each direction, according to the method explained 
in Paper II. 

— First evaluation of level populations, in the optically 
thin approximation. By averaging the incident intensi- 
ties, we obtain mean intensities in the different transi- 
tions of hydrogen and helium. From these intensities and 
physical parameters, we compute the radiative and col- 
lisional transition rates. Then, we solve statistical equi- 
librium equations to obtain atomic level populations at 
each point of the (r, ip) mesh. Since the transition rates 
depend on electron density, it is necessary to iterate. We 
start from an arbitrary value of lo (e. g. lo = 0.5) and, 
using Eqs. (3) and (5), successively deduce: 



A H 



(1 + A Hc + u) k T 



A He = A H A Ho 



and 



N 



A H lo 



(6) 



(7) 



(8) 



After computation of transition rates and solution of 
statistical equilibrium equations, we compute Ahii, 
AhoII, and Ah c iii by adding the populations of indi- 
vidual levels together, and deduce a new value of lo by 



A 



HII 



N\ 



Hell 



2 A, 



HcIII 



A H 



(9) 



These operations are repeated until convergence. 
Full iterations with radiative transfer. This is the main 
part of the computation. The external scheme is similar 
to that of the preceding step, with a variable lo con- 
trolling the convergence of iterations but, in the mean- 
time, the internal intensities for all transitions of hy- 
drogen and helium are recomputed according to the 
principles of NLTE radiative transfer: absorption coeffi- 
cients are derived from atomic level populations (deter- 
mined in the preceding iteration). Then, intensities are 
computed by solving the transfer equation along each 
ray and integrating with respect to direction and fre- 
quency. At the same time, the diagonal terms of the A 
operator are calculated by the method of Rybicki and 
Hummer (1991). The formulae appropriate to the cylin- 
drical geometry are given in Paper II. The new inten- 
sities and the diagonal A coefficients are used to form 
preconditioned statistical equilibrium equations, similar 
to those of Werner and Husfeld (1985). These equations 
are solved to obtain new level populations. Generally, 
one radiative transfer iteration for hydrogen and one 
for helium, between two iterations on lo, are sufficient to 
obtain a good convergence. In a few cases, it was neces- 
sary to perform two radiative transfer iterations for one 
lo iteration. These cases are indicated by an asterisk in 
the last column of Table 1. 
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name 


P 9 (dyn cm~ 2 ) 


T(K) 


A Ha 


vt (km s x ) 


D (km) 


nb. iter. 


tl 


0.1 


6000 


0.1 


5 


1000 


91* 


t2 


0.1 


8000 


0.1 


5 


1000 


81 


t3 


0.1 


10000 


0.1 


5 


1000 


76 


t4 


0.1 


15000 


0.1 


5 


1000 


56* 


t5 


0.1 


20000 


0.1 


5 


1000 


71* 


t6 


0.1 


30000 


0.1 


5 


1000 


43 


t7 


0.1 


40000 


0.1 


5 


1000 


5 


t8 


0.1 


50000 


0.1 


5 


1000 


5 


t9 


0.1 


65000 


0.1 


5 


1000 


5 


tio 


0.1 


80000 


0.1 


5 


1000 


5 


til 


0.1 


100000 


0.1 


5 


1000 


5 


Pi 


0.02 


var. 


0.1 


5 


2000 


45 


P 2 


0.03 


var. 


0.1 


5 


2000 


38 


p3 


0.05 


var. 


0.1 


5 


2000 


33 


p4 (or a4) 


0.1 


var. 


0.1 


5 


2000 


42 


P 5 


0.2 


var. 


0.1 


5 


2000 


40 


P 6 


0.3 


var. 


0.1 


5 


2000 


39 


P7 


0.5 


var. 


0.1 


5 


2000 


51 


al 


0.1 


var. 


0.01 


5 


2000 


41 


a2 


0.1 


var. 


0.02 


5 


2000 


41 


a3 


0.1 


var. 


0.05 


5 


2000 


42 


a5 


0.1 


var. 


0.15 


5 


2000 


42 


a6 


0.1 


var. 


0.20 


5 


2000 


43 


a7 


0.1 


var. 


0.30 


5 


2000 


43 



2.3. Models 

The model cylinders are defined by a series of geometrical 
and physical parameters. The geometrical ones, diameter 
D, altitude H and inclination a, have been defined above. 
Physical parameters include: gas pressure P g , temperature 
T, microturbulent velocity Vt, and relative helium abun- 
dance Aho- The code allows the definition of each of these 
parameters as a function of r and ip, but this possibility is 
not used in practice, except for the temperature. A sum- 
mary of models used in the present paper is displayed in 
Table 1. First, there is a series of isothermal models ("tl" 
to "til") with temperature ranging from 6000 to 10 5 K. 
Their diameter is fixed to 1000 km, their pressure to 0.1 
dyn cm~ 2 , helium abundance to 0.1, and microturbulent 
velocity to 5 km s _1 . Other models have a prescribed tem- 
perature variation T(r), which comes from Paper III, and 
is represented in the upper part of Fig. 1 (models using 
this temperature distribution are indicated by "var." in the 
third column of Table 1). The mean model "p4" (or "a4") 
is similar to that of Paper III. In the series "pi" to "p7", 
we change the pressure while keeping the other parame- 
ters constant. In the series "al" to "a7", we investigate the 
effects of helium abundance. The last column of Table 1 
indicates the number of iterations necessary to achieve the 
convergence on cu, the criterion being fixed to 10~ 6 and the 
minimum iteration number to 5. 

3. Ionization 

The introduction of helium ionization in models allows the 
electron density to be greater than the hydrogen density. 
Let £ = (Nhii/Nh) be the ionization ratio for hydrogen, 
and similarly rjx = (N HeII /N He ) and i] 2 = (N HeIII /N He ) 
be the corresponding ratios for the two stages of helium 
ionization. Equation (9) becomes 

oj = £ + A He 771 + 2 A He t] 2 . (10) 



The ionization ratios £, 771 and 772 are principally controlled 
by temperature. However, near the surface of the structure, 
the influence of incident radiation becomes more and more 
important and tends to moderate the effects of tempera- 
ture. 

3.1. Model with temperature gradient 

In the standard model " p4" , the temperature increases with 
the distance to the axis, as shown in the upper part of Fig. 1. 
In the core of the cylinder at 6000 K, helium is essentially 
neutral and hydrogen is weakly ionized, with a ratio £ be- 
tween 0.2 and 0.3. The first ionization potential of hydrogen 
being lower than that of helium, £ begins to increase first 
when moving towards exterior. This major change in hydro- 
gen ionization occurs between 0.5i? and 0.7i?. For r > 0.7R, 
which corresponds to temperatures T > 20000 K, hydrogen 
is almost completely ionized. The first stage of ionization 
for helium (ratio 7/1) occurs between (r/R) = 0.7 and 0.8, 
i.e. temperatures between 20000 and 35000 K. Finally, the 
second stage of helium ionization starts near (r/R) = 0.9, 
which corresponds to a temperature of 60000 K. These vari- 
ations of £, ?7i and 772 are represented in the middle part 
of Fig. 1. The consequences for electron density are shown 
in the lower part of Fig. 1, with a comparison to the case 
(Paper III) where helium ionization was neglected. The two 
curves practically coincide up to T = 20000 K, since hydro- 
gen is the only electron contributor there. Between 30000 
and 60000 K, oj remains close to 1.1, the main part of he- 
lium atoms being singly ionized. Near the surface, u tends 
to 1.2, that is the possible maximum value since we have 
assumed Ah & =0.1. 

3.2. Isothermal models 

In the standard model, the variations of temperature with r 
masks the effects of incident radiation on ionization. To dis- 
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0.9 1.0 




Fig. 1. Variations of temperature T and population ratios 
with the distance to the axis (r) , for the model " p4" at the 
foot of the loop (a — 0). Abscissae: distance to axis rela- 
tive to the total radius R. Top: temperature. Middle: ion- 
ization ratios for hydrogen (£: dot-dashed line) and helium 
(rji: dashed line; 772: continuous line). Bottom: electron-to- 
hydrogen ratio u> (dashed line: model assuming neutral he- 
lium; continuous line: model with both hydrogen and he- 
lium ionization). 



tinguish between thermal and radiative effects, it is neces- 
sary to examine the variations of ionization within isother- 
mal models. We consider the series of "t" models described 
in Table 1. Here, the cylinders have horizontal axes, but 
the number densities are averaged over ip to obtain a sin- 
gle ionization curve (£(r), r)i(r) or 772(f)) for each model 
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Fig. 2. Variations of the ionization ratios £ (top) and 771 
(bottom) with r for different isothermal models (for clarity, 
all models are not represented) . Number densities are aver- 
aged with respect to the azimuth ip. Symbols for hydrogen: 
open circles: T — 6000 K; full circles: 8000 K; open squares: 
10000 K; full squares: 15000 K; crosses: 20000 K; continu- 
ous line: 50000 K. Same symbols for helium, plus: dotted 
line: 30000 K; dashed line: 65000 K; dot-dashed line: 80000 
K. 



and each ion species. Figure 2 represents the variations of 
ionization of hydrogen and helium for the different mod- 
els of this series. It appears that the curves £(r) and f]i(r) 
are generally flat, which confirms that the effects of tem- 
perature dominate those of incident radiation. This is es- 
pecially true for helium. For hydrogen, the £(r) curves ex- 
hibit some inflexions near the surface. For the cool models 
at 6000 and 8000 K, the incident radiation tends to in- 
crease hydrogen ionization near the surface. In contrast, 
a decrease appears near the surface for models at 10000 
and 15000 K. For temperatures higher than 20000 K, the 
curves £(r) are flat: hydrogen is almost completely ionized, 
so that the cylinders are optically thin in all transitions of 
the hydrogen atom, and consequently insensitive to inci- 
dent radiation. Hydrogen is approximately half-ionized for 
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T = 10000 K. The curves r?i(r) corresponding to the first 
ionization stage of helium are very flat, but a slight decrease 
near the edge may be observed for the model at 20000 K. 
The principal change of ionization for helium occurs be- 
tween 20000 K (jji w 0.3) and 30000 K (771 w 0.9). Between 
30000 and 65000 K, more than 90% of helium atoms are 
in the first stage of ionization (nearly 100% at 50000 K). 
At higher temperatures, rji decreases while 772 increases. At 
80000 K, the numbers of Hell and He III ions are approxi- 
mately equal. At 100000 K, 771 is lower than 0.2, and conse- 
quently 772 greater than 0.8. The electron-to-hydrogen ratio 
u> (not represented) may be easily deduced from the curves 
of Fig. 2, according to Eq. (10). At low temperatures, ui is 
nearly equal to £. In contrast, above 50000 K, there are no 
longer neutral atoms, so that £ = 1 and 772 = 1 — t?i , which 
gives oj w 1 + 0.1 x (2 - 771). 

4. Formation of helium lines 

The computations described in Sect. 2 provide us with 
absorption coefficients and source functions in the different 
lines of hydrogen and helium, which may be subsequently 
used to calculate the intensities emerging from the cylinder. 
The introduction of helium ionization has little influence 
on hydrogen lines, which are formed in regions where 
helium is essentially neutral. Since hydrogen lines have 
been treated in preceding papers, we concentrate here on 
helium lines, and use the standard model "p4" to study 
their formation. 

The absorption coefficients and source functions for 
the different values of r and ip, are computed for 7 values 
of the inclination, from a = to a = 90°. Then, using 
an interpolation procedure described in Paper III, we 
construct a semi-toric loop model and obtain the inten- 
sities emitted towards the observer. Figure 3 shows the 
intensities, integrated over frequency, for a few lines chosen 
among the most important of the helium spectrum. The 
images corresponding to the two optical or infrared lines, 
say 10830 and 5876 A, look very similar to each other. 
The consideration of absorption coefficients indicates that 
the cylinder is optically thin for these two transitions (or 
marginally thick for 10830 at line center). In contrast, 
the two ultraviolet lines, He I 584 and He II 304 A, are 
definitely optically thick. The structure seems broader in 
He II than in He I lines, which is explained below. 

To discuss the process of formation, it is sufficient 
to consider the top of the loop, equivalent to a cylinder 
with horizontal axis. For an observer looking at the cylin- 
der in a direction normal to the axis, the path of photons 
is represented on Fig. 4. The ray enters the cylinder at 
point A and exits at point B. If R is the radius of the 
cylinder and r± the distance between the axis and the ray, 
the abscissae of A and B are — s max and s max , respectively, 
with 



cylinder 
section 



■ _L- 



(11) 



Let k v and S v be the absorption coefficient and the source 
function, respectively, at frequency v. From the transfer 
equation 

dh 

ds 
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Fig. 5. Total emergent intensities computed from model 
"p4" , for different helium lines vs. position across the cylin- 
der. Abscissae are in megameters and grow as the distance 
from the Sun. Ordinates are intensities in erg cm~ 2 s _1 
sr" 1 . Dashed line: He I 10830 A; dotted line: He I 5876 A; 
dot-dashed line: He I 584 A; continuous line: He II 304 A. 



Kv{Si/ Iv) 



(12) 



and the boundary condition 

Iv( Smax) — , 

we deduce the emergent intensity 

/Smax 
«„(«) S v (s) e-^^l ds , 
"Smax 



(13) 



(14) 
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Fig. 3. Emission of the loop model "p4" in several lines of helium: He I 10830 A (top, left); He I 584 A (top, right); He I 
5876 A (bottom, left); He II 304 A (bottom, right). Frequency- integrated intensities are normalized to the maximum 
value of each image. Horizontal and vertical coordinates indicate distances in megamctcrs. 



with the optical thickness between the running point M 
(abscissa s) and B 



r„(a) 



k v (s') ds' 



(15) 



The total emergent intensity for the line under considera- 
tion is then 



-/ 

J line 



^(w) dv 



(16) 



with 

C v {s) = k v (s) S v {s) e 
Thus, Eq. (16) becomes 

7i ine = / C(s) ds 

with 



As) 



These frequency-integrated intensities have been computed 
as functions of r and ip for different helium lines. They are 
shown in Fig. 5. For all these lines, we observe a global 
decrease of intensities from the lower to the upper edge, 
which is due to the decrease of incident radiation. However, 
several differences may be noticed. Concerning ultraviolet 
optically thick lines, the cylinder appears broader in the 
He II resonance line at 304 A than in the corresponding 
line for He I at 584 A. Besides, the two triplet lines of 
He I under consideration have in common a transversal 
variation with three smooth peaks. 

To locate the origin of emissions inside the cylinder, 
we rewrite Eq. (14) as 



CM 



/ ■ 

J line 



C v (s) dv . 



(18) 



(19) 



(20) 



^i/(Smax) 



/Sn 
-Si 



C v (s) ds 



(17) 



C(s) is the contribution function appropriate for the 
frequency-integrated emergent intensity. This function is 
plotted as shades of gray in Fig. 6 for the three lines 
at 10830, 584 and 304 A, respectively. The fi gure for the 
5876 A line is very similar to that of the 10830 A line. The 
existence of three maxima of emission for the intensity at 
10830 A is due to the existence of two zones in the con- 
tribution function: a central patch and a ring. The patch 
corresponds to radiative processes of emission by low tem- 
perature matter, while the ring corresponds to a range of 
temperature where the maximum excitation of He I is oc- 
curing. The radiative processes of emission that occur in the 
central patch include direct scattering of incident radiation 
and the photoionization-recombination process (Hirayama 
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119711 Zirin 11973)) . Since these two processes are dependent 
on incident radiation, the contribution functions in the cen- 
tral zone decreases with height. In contrast, the ring is pro- 
duced by collisional excitation, so that it is practically in- 
dependent of tp. The contribution function for the 584 A 
line is concentrated in a single ring, without central part: 
since the line is optically thick, the front part of the ring 
only is contributing to intensity. The same is true for the 
304 A line but, in this case, the emitting layer is located at 
high temperature, very close to the surface. For this reason, 
the loop looks broader in this transition than in the He I 
lines. 

5. Influence of physical parameters 

5.1. Temperature 

To study the effects of temperature on emitted intensities, 
we use the series of isothermal models ranging from "tl" 
(6000 K) to "til" (10 5 K). The cylinders under consid- 
eration have a horizontal axis (a = 90°) and the emitted 
intensities are averaged along a transversal direction. This 
spatial averaging is convenient for comparison with most 
observations, where the width of the loops corresponds to 
a small number of pixels. Frequency-integrated intensities, 
for the different models and the four principal lines, 
are displayed in Fig. 7. The He I triplet lines at 10830 
and 5876 A have again a similar behavior: the intensities 
decrease from 6000 to 15000 K, rise from 15000 to 30000 K, 
and decrease again at higher temperatures. As mentioned 
in the preceding section, these lines have two sources of 
emission: radiative (dependent on incident radiation) and 
collisional. At low temperatures, scattering is dominant, 
and the decrease of optical thickness results in a decrease 
of scattered radiation. At higher temperatures, collisional 
excitation grows and, despite the decrease of optical thick- 
ness, the emitted intensity rises. For temperatures greater 
than 30000 K, the effect of opacity decrease (principally 
due to ionization, as may be seen in Fig. 2) dominates that 
of excitation increase, so that the intensity decreases again. 
The intensity variation of the resonance line at 584 A is 
similar, but this line is optically thick at low or medium 
temperatures, so that the slope of the curve between 6000 
and 15000 K is not so steep as that of the triplet lines. 
At high temperatures, the decrease of the 584 A line 
intensity is faster than that of the triplet lines, which may 
be due to the decrease of density. According to the analysis 
of neutral helium line formation by Andretta and Jones 
(|1997p . the intensity ratio (triplet/singlet lines) decreases 
with increased density. The behavior of the resonance 
line of He II is quite different and shows a slow and 
steady increase of intensity with temperature until about 
80000 K, followed by a decrease. Most He II lines reach 
their peak of emission between 5 x 10 4 K and 10 5 K (see 
for instance Laming and Feldman l 1 9931 for the 1640 A line). 

The consideration of line profiles (Fig. 8) gives an- 
other insight on the effects of temperature variations. 
The profiles of triplet lines are nearly gaussian, a typical 
feature of optically thin lines, their widths increasing with 
temperature. The 584 A line exhibits a reversed profile 
from 6000 to 30000 K. At 50000 K, the reversal disappears, 
indicating that the line is optically thin. This is due to the 
fast decrease of neutral helium populations above 3 x 10 4 



K (see Fig. 2). At 10 5 K, the central intensity of the 584 
A line is two orders of magnitude lower than at 5 x 10 4 
K. For the He II line at 304 A, the line center intensity is 
roughly constant until 5 x 10 4 K, and the global increase 
of intensity is due to the progressive broadening of the 
line profile. At these temperatures, the main process of 
emission is the scattering of incident radiation, and the 
broadening of the line is due to the joint effects of opacity 
(increase of 771) and temperature (Doppler effect). Above 
5 x 10 4 K, 7/1 begins to decrease, but collisional excitation 
grows. These two competing effects produce a maximum 
of intensity near 8 x 10 4 K. In Fig. 8, one may compare 
the profiles of the 304 A line at 5 x 10 4 and 10 5 K: the 
higher temperature profile is more intense at line center, 
but the intensity decreases more rapidly in the wings, as a 
consequence of the smaller opacity. 

5.2. Pressure 

The effects of gas pressure on emitted intensities are studied 
by means of a series of models "pi" to "p7". These mod- 
els have the same variation of temperature as the standard 
model, but differ from each other by pressures ranging from 
0.02 to 0.5 dyn cm -2 . Frequency-integrated intensities for 
the principal helium lines emitted by these models are dis- 
played in Fig. 9. It appears that all intensities increase with 
pressure, but the slopes of the curves differ from one line to 
another. For neutral helium lines, the variation may be un- 
derstood by considering the two main processes of emission: 
scattering of incident radiation and collisional excitation. 
In the optically thin case, the first process yields intensities 
proportional to atom populations. The second emission pro- 
cess, collisional excitation, which is proportional not only 
to the emitting atom density, but also to the electron den- 
sity, tends to produce a quadratic variation of emission in 
function of pressure. These considerations apply to the two 
He I lines at 10830 and 5876 A. It is visible in Fig. 9 that 
the emitted intensities for these two lines are proportional 
to pressure from 0.02 to about 0.2 dyn cm -2 , and that the 
slope slightly increases at higher pressures, when collisions 
cease to be negligible. The same considerations apply to the 
584 A line, but the change of slope begins at lower pres- 
sures, around 0.1 dyn cm~ 2 . The case of the He II 304 A line 
is more complicated, since the slope of I(P) is smaller than 
linear below 0.1 dyn cm~ 2 and greater than linear at higher 
pressures. Besides, it is visible on profiles (Fig. 8) that this 
line is optically thick and formed in the outermost part of 
the cylinder (Fig. 6). Whatever the pressure, the incident 
ultraviolet radiation from the Sun ionizes helium near the 
surface of the cylinder and creates a zone which scatters the 
304 A radiation. This part of emission due to scattering is 
nearly constant and constitutes the main contribution at 
low pressures (0.02 to 0.05 dyn cm~ 2 ). In contrast, at high 
pressures (0.2 to 0.5 dyn cm -2 ), collisional excitation be- 
comes dominant and results in a quasi-quadratic variation 
of J(P). 

5.3. Helium abundance 

The abundance of helium is one of the most important 
parameters when analysing helium lines (see for instance 
Andretta et al. 120081 and references therein). Some observa- 
tions, like that of filaments by Gilbert et al. (2007| suggest 
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corresponds to a specific line, as indicated inside. 



that very important changes of helium abundance may oc- 
cur between the top and the base of these objects. To eval- 
uate the importance of this parameter on the emission of 
cylinder threads, we use a series of models ("al" to "a7"), 
with the same pressure and temperature variation as the 
standard model, but helium to hydrogen ratios ^4ho vary- 
ing from 0.01 to 0.3. The variations of intensities as func- 
tions of An c arc displayed in Fig. 10. The interpretation 
of these curves is relatively easy. As long as the cylinder 
is optically thin in the considered transition, the intensity 
is proportional to Ah c . This is the case for the 10830 and 
5876 lines until Ah c = 0.05. For higher abundances, the 



slopes of the curves slightly decrease as the line centers be- 
gin to saturate. For optically thick lines at 584 and 304 A, 
the intensity is still a growing function of abundance, but 
the slope is significantly lower than that corresponding to 
proportionality. 

6. Conclusion 

At this stage of our project, it is possible to perform a 
modeling of solar coronal loops including the following in- 
gredients: 

— NLTE radiative transfer in cylindrical geometry, 
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— 2 dimensions (radius and azimuth), 

— statistical equilibrium of atomic level populations, 

— self-consistent treatment of the two principal chemical 
elements: hydrogen and helium, including ionization, 

— pressure equilibrium. 

The numerical code produces in a single run spectral line 
profiles and intensities for the principal transitions of hy- 
drogen and helium, emitted in all directions. Some develop- 
ments remain to be done, as the implementation of motions 
in helium routines (as was done for hydrogen in Paper V), 
or the treatment of Lyman lines in partial redistribution. 
We also are willing to apply our code to the modeling of 
complex loop systems observed in the solar atmosphere. In 
this scope, we envisage to extend the code to some minor 
chemical species, in function of the contents of observations. 

Acknowledgements. We wish to thank Pctr Hcinzel for useful sugges- 
tions concerning this manuscript. 
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Fig. 6. Contribution functions corresponding to the stan- 
dard model with horizontal axis, for three lines of helium. 
Top: He I 10830 A; Middle: He I 584 A; Bottom: He II 304 
A. Bright zones correspond to the maximum of the func- 
tion. Directions are the same as in Fig. 4 (observer at right, 
Sun below). 
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Fig. 7. Frequency-integrated intensities, averaged over 
position, emitted by isothermal models "tl" to "til". 
Intensities (erg cm -2 s _1 sr _1 ) arc plotted as functions of 
the temperature (K) of the model, in 4 transitions: He I 
10830 A (open circles); He I 584 A (full circles); He I 5876 
A (open squares); He II 304 A (full squares). 
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Fig. 10. Frequency-integrated intensities, averaged over 
position, emitted by models "al" to "a7". Abscissae: log 
(abundance ratio An e ). Ordinatcs: intensities in erg cm~ 2 
s _1 sr _1 . Spectral lines: He I 10830 A (open circles); He I 
584 A (full circles); He I 5876 A (open squares); He II 304 
A (full squares). 
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Fig. 9. Frequency-integrated intensities, averaged over po- 
sition, emitted by models "pi" to "p7" . Abscissae: gas pres- 
sure inside the cylinder (dyn cm -2 ). Ordinates: intensities 
in erg cm -2 s _1 sr _1 . Spectral lines: He 1 10830 A (open cir- 
cles); He I 584 A (full circles); He I 5876 A (open squares); 
He II 304 A (full squares). 



